Astronomy & Astrophysics manuscript no. 5768 


© ESO 2008 


February 5, 2008 





On the asymptotic giant brancli star origin 
of peculiar spinel grain 0C2 



\o ■ 

O . M. Lugaro\ A. I. Karakas^, L. R. Nittler^, 

O ■ C. M. O'D. Alexander^ P. Hoppe"^, C. Iliadis^ and J. C. Lattanzio^ 

CN ■ 



O ! ' Sterrenkundig Instituut, University of Utrecht, Postbus 80000, 3508 TA Utrecht, The Netherlands 
' ^ Origins Institute, Department of Physics & Astronomy, McMaster University, Hamilton ON, Canada 
Carnegie Institution of Washington, Department of Terrestrial Magnetism Washington DC 20015, USA 
Max-Planck-Institute for Chemistry, P.O. Box 3060, 55020 Mainz, Germany 
' ' ■ ^ Department of Physics and Astronomy, 176 Phillips Hall, University of North Carolina, Chapel Hill, NC 27599-3255, USA 
I _ * Centre for Stellar and Planetary Astrophysics, School of Mathematical Sciences, Monash University, Victoria 3800, Australia 

> : 

' Received June -, 2006; accepted — , 2006 

\o : 

■ ABSTRACT 

o ■ 



o 



Context. Microscopic presolar grains extracted from primitive meteorites have extremely anomalous isotopic compositions revealing the stellar 
origin of these grains. Multiple elements in single presolar grains can be analysed with sensitive mass spectrometers, providing precise sets of 
isotopic compositions to be matched by theoretical models of stellar evolution and nucleosynthesis. 
' Aims. The composition of presolar spinel grain OC2 is different from that of all other presolar spinel grains. In particular, large excesses of the 
Oh' heavy Mg isotopes are present and thus an origin from an intermediate-mass (IM) asymptotic giant branch (AGB) star was previously proposed 
^ , for this grain. We discuss the O, Mg, Al, Cr and Fe isotopic compositions of presolar spinel grain 0C2 and compare them to theoretical 
5_j ' predictions. 

Methods. We use detailed models of the evolution and nucleosynthesis of AGB stars of different masses and metallicities to compare to the 
d : composition of grain OC2. We analyse the uncertainties related to nuclear reaction rates and also discuss stellar model uncertainties. 
L| ' Results. We show that the isotopic composition of O, Mg and Al in 0C2 could be the signature of an AGB star of IM and metallicity close to 
. ^ , solar experiencing hot bottom burning, or of an AGB star of low mass (LM) and low metallicity (=^0.004) suffering very efficient cool bottom 
■ processing. Large measurement uncertainty in the Fe isotopic composition prevents us from discriminating which model better represents the 
?H parent star of OC2. However, the Cr isotopic composition of the grain favors an origin in an IM-AGB star of metallicity close to solar. 

Conclusions. Our IM-AGB models produce a self-consistent solution to match the composition of 0C2 within the uncertainties related to 
reaction rates. Within this solution we predict that the '*0(p, 7)'^F and the "0(p, q')'"'N reaction rates should be close to their lower and upper 
limits, respectively. By finding more grains like OC2 and by precisely measuring their Fe and Cr isotopic compositions, it may be possible in 
the future to derive constraints on massive AGB models from the study of presolar grains. 

Key words. Nuclear reactions, nucleosynthesis, abundances - Meteors, meteoroids - Stars: AGB and post- AGB 



1. Introduction 

Presolar grains were born in circumstellar regions around an- 
cient stars, ejected into the interstellar medium, preserved dur- 
ing the formation of the Solar System, and trapped inside primi- 
tive meteorites from which they are now extracted and analysed 
by various microanalytical techniques. Their isotopic com- 
positions are extremely anomalous compared to those found 
in the bulk of materials formed in the Solar System (e.g.. 
Earth, meteorites, etc.). They represent a detailed record of the 
compositions of their parent stars, and, as such, provide ma- 
jor constraints for models of stellar structure and nu cleosyn- 
thesis and of the chemical evolution of the Galaxy dZinneil 
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Il998t IClavton & Nitdeil Eooi iLugarol l2005l) . Until recently, 
Mg-rich presolar grains have been thought to be rare, as the 
dominant identified carbide (SiC, graphite) and oxide (corun- 
dum: AI2O3) phases have low Mg contents. The situation 
changed with the advent of a new type of ion microprobe 
with improved spatial resolut i on an d sensitivity, the Cameca 
NanoSIMS 50 jHiUion et al.L Il994 . Using the NanoSIMS, 
IZinner et alJ ( l2003h found that presolar spinel is, in fact, more 
abundant in meteorites than is presolar corundum, but has a 
finer grain size, with most grains having diameters less than 1 
jum. 

A presentation and discussion of the properties and the 
compositions of the more than 3 00 presolar spinel g rains 
reported so far can be found in IZinner et alJ ( l2005h and 
iNguven. Zinner & LewisI (l2003h . Most presolar oxide as well 
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Fig. 1. Scanning electron microscope image of presolar spinel 
grain 0C2. This 800nm grain is sitting on a gold pedestal, fol- 
lowing the ion probe isotopic analysis, because the gold sub- 
strate sputters faster than the grain does. 
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Fig. 2. Mg isotopes in presolar spinel and corundum grains 
( |Nittler_et al . 1997; Choi et al , 1999; Zinn er et aL. 20ol . 
Dashed lines indicate solar isotopic ratios; errors are l-cr. 
Spinel grain OC2 has larger excesses of ^^Mg and ^^Mg than 
seen in other presolar spinel grains. Measured 6(-^Mg/-'^Mg) 
values for corundum extend up to 1.7x10*', due to high initial 
contents of -^Al and high Al/Mg ratios, but no corundum grain 
has a ^^Mg enrichment similar to that of OC2. 



as silicate grains show enhancements in '^O and depletions in 
'^O similarly to what is observed in red giant and low-mass 
(LM) AGB stars, p ointing to such an origin for the major- 
ity of these grains jNittler et all Il997l) . These compositions 
are explained by the effect of the first dredge-up and extra- 
mixing processes. We focus in this paper on the composition 
of a single extraordinary presolar spinel grain, named OC2. 
OC2, of size a;0.8 fim (Fig. [!), was identified in a mixed 
acid residue of the Semarkona, Krymka and Bishunpur un- 
equilibrated ordina ry chondrites jNittl er & Alexande 3 Il999t 
IZinner et al.l l2005l) . The isotopic composition of this grain 
(TableQ) is peculiar among presolar oxide grains. Its most re- 
markable feature are large excesses of the heavy Mg isotopes: 
^^Mg and ^^Mg are enriched with respect to solar composition 
by 43% and 1 17%, respectively (Fig.|2}. The Mg isotopic com- 
position is represented in Table^by the 6 notation; 



5('Mg/24Mg) 



('Mg/24Mg), 



'measured 



('Mg/24Mg) 



solar 



1 X 1000, 



which we will use throughout this paper, where the abundant 
■^'^Mg, which represents about 79% of all magnesium in the 
Solar System, is used as the reference isotope. 

The O isotopic composition of grain OC2 (Fig. |3} is also 
quite extreme with respect to typical compositions of presolar 
oxide grains. While its '^0/'*0 ratio is enriched by a factor 
of 3.3 relative to solar, its "^O/'^O ratio is at least 26 times 
smaller than that of the Solar System. Presolar oxides with 
lar ge '^O deple t ions have be en classified as "Group 2" grains 
bv iNittler et alJ ill 994 Il997l) . but these are in fact quite rare. 
Of some 600 known presolar corundum, hibonite (CaAliaOig) 




18Q/16Q 

Fig. 3. O isotopic ratios measured in presolar corundum, hi 
bonite, spinel, and silicate grains ([Nittle^ Il99 
1997 | IChoi etall Il998t Messenger et all 120031 



Nittler et al 



Zinner et al 



2003; Mostefaoui & Hoppel l2004t iNguven & Zinneit 12004 : 
Nittler et alj.i2005i) . Grain OC2 is a member of the rare Group 2 



class of highly ^O-depleted grains. 



and spinel grains, only 10 have '^O/'^O ratios lower than 10 
similar to that of OC2 (Fig.g). 

The Cr isotopic ratios were also measured in grain OC2, 
as well as the ^^Fe/^^Fe ratio. From Table ^we see that these 
values are typically solar within the error bars. 

The larg e ^^Mg and ^*Mg excesses observed in OC2 led 
IZinner et'all (l2005h to propose for this grain an origin in an 
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Table 1. The measured composition of presolar spinel grain 
OC2. 





OC2 


solar ratio 


<5(--'Mg/-*Mg) 


433.0 ± 10.0 


0.1266 


<5(-^Mg/«Mg) 


1170.0 ± 15.0 


0.1393 


Al/Mg 


2.18 ±0.39 


0.079 


17Q/16Q 


1.25 ± 0.07 xlO^' 


3.83 xlO-"* 


18Q/160 


6.94 ± 1.34x10-5" 


2.00 xlO-3 


5("Fe/5«Fe) 


170 ± 191 


0.0231 


<5(^"Cr/52Cr) 


26 ± 71 


0.0519 


(5(53Cr/52Cr) 


-56 ± 45 


0.1135 


5(^^Cr/^'Cr) 


102 ± 117 


0.0282 



All errors are given at Icr. 

"Quoted here is the experimental statistical uncertainty, which does 
not take into account the possible effect of pollution of terrestrial ma- 
terial discussed in Sec.|3| 

intermediate-mass (IM)-AGB star, wit h mass betwee n 4-7 
Mq and underging hot bottom burning dHerwiglEooi . instead 
of the usual low-mass origin assigned to the majority of the 
other presolar oxide grains. As for the late evolutionary phases 
of most stellar types, strong mass-loss and dust formation are 
observed around massive AGB stars. In fact, IM-AGB stars are 
part of the family of OH/IR stars, i.e. O-rich AGB stars ex- 
tremely bright in the infrared, and astronomical observations 
indicates that OH/IR s tars are the second most irnportan t source 
of dust in the Galaxy (I AlexandeJ fl997l: IWhittetl 1 1 992h . 

During the AGB phase, the H and He-burning shells are 
activated alternately in the deep layers of the star while the 
extended H-rich envelope loses mass through strong stellar 
winds. During each He-burning episode, the sudden thermal 
runaway triggers the whole region between the H- and the 
He-shell (He intershell) to become convective. After a thermal 
pulse (TP) quenches, the base of the convective envelope can 
sink into the He intershell (third dredge-up, TDU), thus car- 
rying to the surface nuclear processed material (see iHerwid 
I2OO5I for a review on AGB stars). 

The nucleosynthesis occurring in IM-AGB stars is different 
from that occurring in LM-AGB stars because the temperature 
during thermal pulses can exceed ^ 350 million degrees and 
thus the heavy Mg isotopes can be produced by a-capture re- 
actions on -^Ne, which is abundantly present in the He inter- 
shell due to the conversion of ''^N from H-burning ashes into 
^-Ne by d ouble a captures durin g the early phases of a ther- 
mal pulse jKarakas et aill2006bl) . The heavy Mg isotopes are 
then mixed to the envelope by the following TDU episode. In 
IM-AGB stars proton captures also occur at the base of the con- 
vective envelope (hot bottom burning, HBB). The MgAl chain 
is activated resulting typically in the destruction of ^^Mg and 
the production of ^^Al and ^^Mg. If the temperature exceeds 
^ 80 million degrees ^"*Mg also suffers proton captures (see 
detailed discussion an d models bv lKarakas & Lattanzioll2003t 
iKarakas et al.[l2006bl) . HBB also greatly affects the O isotopic 
composition, resulting in the production of '^O and the destruc- 
tion of '^O, as observed in OC2. HBB can als o prevent the sur- 
face o f the star from becoming carbon-rich jBoothrovd et al.L 
Il993h . by converting the dredged-up '^C to '"^N. Thus, IM- 



AGB stars might be expected to preferentially form oxide 
rather than carbonaceous phases like SiC and graphite, which 
require C>0. However, the formation of spinel cannot be com- 
pletely ru led out for C/O rati os slightly above unity (see dis- 
cussion in IZinneret al.Ll2005l) . 

Stellar sources of dust in the Galaxy other than AGB stars 
(see, e.g.. Table 1 of lAlexandej 1 19971) seem unlikely to be 
the site of origin of grain OC2. The ejecta of supernova ex- 
plosions are enriched in ^^Mg, ^^Mg and ^^Al when consider- 
ing solar metallicity models. However, contrarily to the com- 
position of OC2, '^ O and '^O are produced, while '^O is typi- 
cally much depleted jRauscher et al . l'2002^'Limongi & Chieiffil 
Eo03). Wolf-Rayet stars, of whi ch the WC typ e are a lso ob- 
served to generate dust (see e.g Iwilhams et al.t 1 19871) . could 
also produce the O isotopic composition of grain OC2 together 
with excesses in "^ Al at a time just befo re the transition from 
WN to WC occurs jArnouldet al.[ll997l) . However, no produc- 
tion of ^^Mg is expected, at least until core He burning starts 
and the star moves into the WC/WO phases, at which point also 
large depletions of '^O and possible enhancements in '^O are 
predicted. Nova nucleosynthesis is predicted to produce high 
enhancements in ^^Mg, '^Mg and ^^Al, however, these are typ- 
icall y accompanied by high enhancements al so in both ' ^O an d 
1^0 llose et aLlEool . Only model COl of I.Tose et alJ t004 . 
a CO nova of low mass 0.6 Mq, produces a deficit in '^O com- 
parable to that observed in grain OC2. However, the '^O/'^O 
in this model is three times higher than in OC2, and the ex- 
cesses in the heavy Mg isotopes are: 5(^^Mg/^''Mg)=836 and 
5(^''Mg/^'*Mg)= 19062 (without adding up the contribution of 
-*'A1, with ^''A1/2^A1=0.006)! Even if mixing between this type 
of nova material and material of solar-like composition were 
invoked, it would be impossible to avoid producing ^^Mg ex- 
cesses much higher than those observed in grain OC2 at the 
time when the ^^Mg/^'^Mg ratio is matched. 

The aim of this paper is to analyse the O, Mg, Al, Cr and 
Fe isotopic compositions predicted by detailed models of AGB 
stars of different masses and metallicities and discuss them in 
the light of the precise measurements of the composition of 
grain OC2. In this way, we can test the idea that grain OC2 
originated in an IM-AGB star. There are few direct constraints 
available to test theoretical models of these stars mostly owing 
to the fact that they are much rarer than their lower-mass coun- 
terparts. There are observations available for Li, C, O, '^C/'^C 
and heavy elements abundances, in particular for AGB stars in 
the M agellanic Cl ouds (see e.g. IWoodetallll9 83 HPle7 etal.L 
dH [Smith et all ll995[lvan Loon et allll999l) . The Li, C and 
O abundances and the '^C/'"^C ratios have been used as tests 
for the occurrence of HBB. Stellar models are able to ex- 
plain the fact that the majority of AGB stars of high lumi- 
nosity are O rich, as well as to reproduce the observed lo\y 
'^C/'^C ratios and high Li a bundances (e.g. iBoothrovd et all 
ll993HMazzitelli et al.Lll999l) . 

The paper is structured as follows: in Sect. |2]we describe 
the numerical method and the structure features of the stellar 
models presented in this paper In Sect.|3]we compare the com- 
position of grain 0C2 to the compositions derived from our 
stellar models and discuss the model uncertainties. In Sect. 0] 
we outline our conclusions. 
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2. Methods and models 

The evolution of stars of different masses and metallicities was 
computed from the zero-age main sequence to near the tip of 
the AGB using the Monash version of the Mt . Stromlo Stellar 
Structure code (see iFrost & Lattanziot Il996l and references 
therein for details). The models of IM-AGB stars that are used 
throughout this paper have masses of 5 Mq with Z = 0.02 and 
0.008 and 6.5 Mq with Z = 0.02 and 0.012. The Z = 0.02 
com putations were assumed to have initial abundances taken 
from lAnders & Gr evesse' ('19891, where as the Z - 0.01 2 mod- 
els have initial abundances ta ken from lAsplund et alJ (l2005l) 
for elemental abundances and iLoddersI ( l2003l) for initial iso- 
topic ratios. Most presolar grains recovered from meteorites 
have come from LM-AGB stars so we also present results from 
one low-mass low-Z AGB model with 2.5 Msun and Z = 0.004. 
Models of AGB stars with Z < 0.004 (or Z = 0.004 in the case 
of IM-AGB stars) are not presented because presolar grains 
most likely originated in field stars from the solar neighbor- 
hood. While we do not have any reliable observational infor- 
mation on inter stellar dust lifetime s, theoretical estimates give 
values < 1 Gyr Clones et allll997l) . 

All models include mass loss on t he AG B according to the 
prescription of IVassiUadis & WoodI lll993l) . show deep TDU 
and with the exception of the 2.5 M© case, HBB, which pre- 
vents them from becoming C rich. For the 5 Mq Z = 0.02 
model we tested the effect of using the mass-loss rate from 
Ivan Loon et al.l ( l2005l) . We find that the two formulas produce 
the same mass-loss rate in the final phases of the evolution (the 
last 8 o ut of 24 pulses), which is consistent given that the for- 
mula of lvan Loon et al.l (l2005l) is valid for the superwind phase 
of extreme mass loss in O-rich AGB stars. In Table|2we present 
some features of the stellar models, where M5Z02 denotes the 
5 Mq model with Z = 0.02. We include the number of TPs, the 
core mass at the first TP, the total amount of matter dredged- 
up into the envelope, an d the maximum dred ge-up effi- 

ciency, A (for the definition, see lKarakas et alll2002l) . where all 
masses are given in solar units. We also present the maximum 
temperature at the base of the convective envelope, T^^^, and in 
the He-intershell, T^™, both in units of 10^ K; and the core and 
envelope mass at the final time-step. In all IM-AGB models, 
HBB was shut off as a consequence of mass loss. Once the en- 
velope mass drops below about 1 .5 Mq the temperature rapidly 
drops below that re quired to sustain proton capture nucleosyn- 
thesis jKarakas et al.n2006b.) . Dredge-up may continue, adding 
further '^C to the envelope an d, depending on the initial metal- 
licity, i ncrease th e C/O to > 1 ('prost e t al.[ll998HKarakas et all 
l2006bll . From iKarakas et al J (.2006b.) . however, this only occurs 
in IM-AGB stars with Z < 0.004. 

We performed detailed nucleosynthesis calculations using a 
post-processing code which includes 74 species, 506 reactions 
and time-depen dent diffusive mixing in all convective zones 
llCannonlll99i . Det ails on the po s t-proce ssing code and net- 
work ar e outlined in iLugaro et al.l (l2004bl) and iKarakas et al.l 
ll2006hll and will not be r epeated here. Details on the reac- 
tion rates can be found in iLugaro et al ] (l2004bl) . where it is 
outlined which of the proton, a and neutron capture reac- 
tion rates we have updated (from the REACLIB data tables 



iThielemann et al.|. Il986l 1991 version), according to the lat- 
est experimental r esults. More recent up dates include the ^^Ne 
+a reaction rates (^ Karakas et al.[E()06bl) and the NACRE rates 
(|Angulo et al., 1999| NACRE) for the proton captures rates of 
the NeNa and MgAl chains. The case of the '^0(/?, Qr)'''N reac- 
tion, of particular importance for the present work, is discussed 
in Sect.O 

For the 5 Mq Z = 0.008 and 2.5 Mq Z = 0.004 models, 
we assume an initial (^-enhanced mixture for the elements O , 
Ne, Mg, Si and S typical of thin-disk stars jReddv et al.Ll2003h . 
The initial composition of the O and Mg isotopes are deter- 
mined by the o'-enhancement of '^O and ~^Mg, so that [O/Fe] 
= +0.4 and [Mg/Fe] = H-0.27 at [Fe/H] = -1, while the abun- 
dances of the neutron-rich isotopes are scaled with the initial 
metallicity. Thus, the '^O/'^O and '**0/'^0 ratios are lower than 
their solar values (see Table [0, by factors 1.44 and 1.91 at 
Z = 0.008 and 0.004, respectively. With this choice we ob- 
tai n, at any given metallicity, higher O ratios than suggested 
by iTimmes et al ] 11295) by following the Galactic chemical 
evolution, and where the O ratios are scaled with the metal- 
licity. However, there are still large uncertaintie s in the evolu- 
tion of the O isotopic ratios with time ( Prantzo s et al.l 11994 
iRomano & MatteuccilbOOsh . Most importantly for the discus- 
sion here, the effect of the first and second dredge-up mixing 
events and HBB almost completely erases any record of the 
initial O ratios at the surface during the TP- AGB phase. 

3. Comparison of model predictions and the 
composition of 0C2 

In Fig. 0] the Mg isotopic composition of grain OC2 is com- 
pared to that predicted by our calculations of IM-AGB models. 
The predictions shown in this plot are calculated by including 
the abundance of -*A1 multiplied by 25, taking into account the 
fact that Al was incorporated in spinel grains during their for- 
mation approximately 25 times more preferentially than Mg, 
given that stoichiometric spinel by definition is MgAl204, i.e. 
it has Al/Mg^2, wh ile this ratio is 0.079 in the Solar System 
dZinner et alll2005l) . We have explored the uncertainty related 
to this ^^Al multiplication factor by considering the Al/Mg ra- 
tios predicted at the stellar surface in our models at the time 
when 5(-^Mg/-'*Mg)oc2 is matched, which are 5% to 20% 
lower than in the Solar System, and the error bars of the Al/Mg 
ratio measured in OC2 (see Table [Q. The resulting uncertain- 
ties are represented by the error bars connected to the predic- 
tion lines in Fig.|3] 

First, let us consider the 5 M© Z = 0.02 model, which does 
not reach the excess in ^^Mg shown by OC2. The ^^Mg/-^Mg 
ratio is determined by the combined effect of the He intershell 
temperatures and the total amount of material mixed to the sur- 
face by the TDU. An increase in either of those two quantities 
would lead to higher ^^Mg/^'*Mg ratios, so the OC2 data indi- 
cate that an IM-AGB parent star of this grain must have expe- 
rienced either temperatures higher than 352 million degrees, or 
a TDU mass higher than 0.05 M©, or both. Second, it appears 
that a narrow range of HBB temperatures is required in order 
to produce enough ^^Al to match the excess at mass 26 shown 
by OC2 at the given 5(^^Mg/^'^Mg) value, avoiding an increase 
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Models: 


M5Z02 


M6.5Z02 


M6.5Z012 


M5Z008 


M2.5Z004 


No. TP 

M,(l) 

'"dred 

^BCE 
ymax 
-"He 
M' 

-'"core 


24 
0.861 
5.027(-2) 
0.961 
64 
352 
1.499 
0.874 


40 
0.951 
4.696(-2) 
0.910 
87 
372 
1.507 
0.963 


51 
0.956 
6.483(-2) 
0.940 
90 
370 
1.389 
0.967 


59 
0.870 
1.745(-1) 
0.952 
81 
366 
1.387 
0.886 


34 
0.602 
1.869(-1) 
0.820 

308 
0.685 
0.673 



o 




j: I I I I I I I I I I I I I I u 

500 1000 1500 

delta (2^Mg + 25x26A|)/24Mg 



Fig. 4. The Mg isotopic compositions of grain OC2 is com- 
pared to our models of IM-AGB stars. The 2cr uncertainties for 
OC2 are roughly within the symbol. Each symbol for model 
predictions represents the composition after a TDU episode. As 
indicated in the x-label the <5(^''Mg/^'*Mg) measured in OC2 is 
compared to prediction lines calculated by including the abun- 
dance of -*A1 multiplied by a factor of 25 to take into account 
the fact that Al is preferentially included in spinel by such fac- 
tor (see text). The uncertainty range in the predictions derived 
from the calculation of this factor are represented by the error 
bars connected to each prediction line at 5(^^Mg/^"'Mg)=433. 
Solar composition is represented by ticked axis at 5 = 0. 

of the i5(^^Mg/^^Mg) value itself. This occurs in the models of 
5 Mq and Z = 0.008 metallicity and 6.5 Mq and Z = 0.02 
metallicity (in this latter case within reaction rate uncertainties, 
see Sect. 13. 1> . where the maximum temperature at the base of 
the convective envelope reaches 8 1 and 87 million degrees, re- 
spectively. For the 5 Mq Z = 0.02 case the temperature is too 
low, < 64 million degrees, to produce ^^Al, while in the 6.5 Mq 
Z = 0.012 case the temperature is too high, < 90 million de- 
grees, so that also ^^Mg is produced by HBB. In principle, it 
is possible that for temperatures somewhat higher than 90 mil- 
lion degrees one could achieve again the composition of OC2 
via HBB, as more ^^Mg is converted into ~^A1. However, these 
temperatures are not achieved in our models, and, moreover, 
they would pose a problem to match the '^O/'^O ratio of the 
grain (see discussion below and in Sect. 13. 1> . 

Note that the 6(-^Mg/-'^Mg) and 6(-^Mg/-^Mg) values for 
the 5 Mq Z - 0.008 model reach values higher than 2000 and 



5000, respectively, outside the range shown in Fig.|3 (For the 
6.5 Mq Z - 0.012 model, also outside the range of the figure, 
5- values reach ^2000 and ^2500, respectively). This could rep- 
resent a difficulty for matching the composition of grain OC2 
using the 5 Mq Z = 0.008 model. It is much more likely for 
a grain to be produced with a composition that corresponds 
to abundances in the envelope at the end of the AGB. This 
is because radio and infrared observations confirm that mass- 
loss rates of AGB stars increase with time and that these stars 
expel a large part of their envelope towards the end of their 
evol ution via a strong superwind, with M ~ few lO ^^ Mo/yr 
(e.g. IVassihadis & WoodI Il993[ Ivan Loon et all l2005l) . lead- 
ing to the formation of a planetary nebula. In the case of our 
5 MqZ - 0.008 model, about 40% of the mass is lost in the 
last few TPs when the Mg isotopic ratios are well above that 
observed in OC2. Instead, the mass lost in between the TDU 
episodes that cover the composition of OC2 is only ^ 0.005 
Mq, hence, the probability of a grain forming with such compo- 
sition is much smaller, ^ 0.1%, even if this occurrence cannot 
be ruled out. From this point of view, the 6.5 Mq Z — 0.02 
model could be favored for the parent star of OC2, since it 
reaches the observed (5(^^Mg/-'*Mg) towards the end of its evo- 
lution and the 5(^*Mg/^'*Mg) value is within reaction rate un- 
certainties (see Sect. 13. U . However, it should be also kept in 
mind that the mass-loss prescription is one of the largest uncer- 
tainties in AGB models, as discussed in Sect. 13.21 

In Fig.|5] the O isotopic composition of grain OC2 is com- 
pared to that predicted by our IM-AGB models, and in Fig.|6] 
the '^O/'^O ratio is plotted as function of the (J(^^Mg/^'*Mg) 
value. These figures show that the '^O/'^O ratio of grain 0C2 
is not matched by any of the r nodels, as they alway s produce too 
high a ratio. As discussed bv lNolIett et"ai] ll2003l). the '^O/'^O 
equilibrium ratio during proton captures is mainly determined 
by the ratio of the rates of the nuclear reactions that produce 
and destroy '^O, i.e. the "'0(p,r)'^F and the '^0(p, q;)''*N re- 
actions, respectively. This ratio reaches a minimum of ^ 0.001 1 
around 50 million degrees and then increases again for higher 
te mperatures reaching 0.008 at 100 million degrees (see Fig. 8 
of lNolIettetaill2003l) . During HBB, the '^O/'*"© equilibrium 
ratio is reached at the base of the convective envelope and the 
envelope material is efficiently replaced by material with the 
'^0/"'0 equihbrium ratio after ~ 5 TPs. The 5 Mq Z = 0.02 
model comes closest to producing the needed ratio because of 
its lower temperature, however, this model is the furthest from 
matching the measured 5(-^Mg/-'*Mg) value. The high temper- 
atures at which HBB operates in the 5 Mq Z - 0.008 and in the 
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Fig. 6. The '^O/'^O ratio is plotted as function of the 
5(25Mg/24Mg) for our IM-AGB models. 

6.5 Mq Z = 0.02 models lead to an even worse match with the 
O isotopic composition of grain OC2. The 6.5 M© Z = 0.012 
model follows the expected trend as the '^0/'*0 is the highest 
in this case. 

As for the '^O/'^O ratio, because of HBB, all the models 
reach ratios of the order of 10"^ - 10"^, much lower than that 
shown by grain OC2, after only about 5 to 15 thermal pulses, 
much earlier than when the required Mg isotopic composition 
is reached. So, it is not possible to match the observed value by 
any of our models, but we note that there is always surface con- 
tamination on sample mounts and residual oxygen in the ion 
microprobe vacuum system. The very low measured "*0/'^0 
ratio for grain OC2 was based on 35 actual counted '^O atoms. 



If the grain actually had an '^O/'^O ratio of zero, this low mea- 
sured '^O signal would correspond to a 2% level of terrestrial 
contamination, which is perfectly reasonable. Thus, we con- 
sider it likely that the true "^O/'^O ratio of OC2 was indeed 
lower, and do not consider the mis-match with HBB models to 
be a major problem. 

3.1. The effect of varying the reaction rates 

As discussed above, the '^0/'*0 ratio only depends on the ra- 
tio of the '^0(/7,7)'^F and the '^©(p, q')''*N reaction rates, so 
it is important to carefully evaluate the uncertainties related to 
these reaction rates in the temperature range of HBB, ^ 60 to 
100 million degrees. The '^0(/9, a)'"'N reaction rate in this tem- 
perature range is almost completely determined by a resonance 
at 70 keV. In the models presented above we h ave use d the rate 
fromiBlackmon et al. ( 1995) and Landre et aLl (ll990l) . which is 
the same as in the NACRE compilation. We have also calcu- 
lated a new rate for this reaction on the basis of the l atest avail- 
able e xperimental information dChafa et al.L l2005t IFox et all 
I2OO5I) . We have found a rate close to NACRE (within a few per- 
cent, for the temperature of interest in this study) with an uncer- 
tainty range of ^ -1-25% and -30%, compared to the NACRE 
ranges of ^ + 33% and -22%. The '^0(/7,7)'^F rate is mostly 
determined by direct capture and the astrophysical S'-factor at 
zero energy is extrap olated from the ava ilable data using a 
given potential model jAngulo et allll99'9l) . At low energy the 
experimental data show large error bars and the NACRE com- 
pilation adopts a 30% uncertainty for the derived S -factor This 
results in lower and upper limits for the rate a factor 1 .43 lower 
and 1.30 higher, respectively, than the recommended value at 
the temperature of interest. As shown in Fig.0 the use of the 
upper limit for the '^0(/?,y)'^F rate together with the lower 
limit for the ^^0(p,ay'^0 rate is predicted to provide a good 
match to OC2, using the same models (5 Mq, Z - 0.008 and 
6.5 Mq, Z - 0.02) that match its Mg isotopic composition. 

The Mg isotopic composition is affected by the uncertain- 
ties in the -^Ne-Ha reaction rates as well as proton-capture reac- 
tion rates. The uncertainties of the ^^Ne-i-g reaction r ates have 
been recently re-evaluated by iKarak as et alJ (l2006bl) and are 
relatively small (of the order of ^60%), while the uncertain- 
ties on the ^^Mg(/7, 7)^''A1 and ^^Al{p, 7)^'Si reaction rates are 
large (a factor of 3 and a factor of 1200, respectively) and af- 
fect, in particular , the production of ^ ^Al. If we apply the uncer- 
tainties found bv llzzardetal]ll2006l) to the ^^Al surface abun- 
dances of our models we find that the lower limit for the -^Al 
yield, is incompatible with the composition of OC2, while the 
upper limit for the ^*A1 yield would allow up to ^ 50% more 
^*A1 production. Combining this effect with the uncertainties 
of the ^^Ne+a reaction rates (i.e. using the upper limit for the 
^^Ne(Q', 7)^^Mg and the lower hmit for the ^^Ne(Q', «)^^Mg re- 
action rates, respectively), we find that our 6.5 Mq, Z = 0.02 
model can also provides a match to the Mg isotopic composi- 
tion of OC2. 

As a final remark we note that our stellar structure cal- 
culations have been performed using a '^N(p,y)'^0 reac- 
tion rate very similar to NACRE. However, direct experimen- 
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Fig. 7. Selected IM-AGB model predictions from Fig. |6] are 
shown, with extra model predictions calculated using the lower 
limit (LL) and the upper limit (UL) for the '*0(/7, y) '^F and the 
^^0{p, a)'*N reaction rates, respectively. 



tal data are now available for this rate down to the t ypical 
H-burning temperatures of AGB stars dRunkle et al.L l2005t 
ILUNA coUaboratiorl Eoo3) . They give a rate 40% lower than 
NACRE. This has been shown to have an irnpact on the struc- 
ture of AGB stars, in particular lWeiss et al.l ll2005l) have shown 
that for a 5 Mq star of solar metallicity the peak luminosity 
in the rmal pul s es is higher, a nd the interpulse duration longer, 
while iHerwi g| ll2004cll2006l) have shown that for a 2 Mq star 
of Z =0.01 the choice of the new ''*N(/7, y)'^0 reaction rate 
results in a more efficient third dredge-up. Further analysis is 
needed to test the effect of this updated reaction rate on nucle- 
osynthesis in IM-AGB stars. 

3.2. The effect of model uncertainties 

One of the main uncertainties in AGB models is the choice 
of mass loss: a definitive a greement on the description of this 
phenomena is still missing' Ve ntura & D'Antonal (12005 al) cal- 
culated models with different choices of mass loss and showed 
that a stronger mass loss leads to less efficient HBB nucleosyn- 
thesis. In this case it may be more difficult to produce enough 
26 Al to match the ^(-''Mg/^^Mg) of OC2. 

Another uncertainty in AGB star models relates to the pos- 
sible mixing of protons from the H-rich envelope into a tiny 
region at the top of He intershell (partial mixing zone, PMZ) 
at the end of each TDU episode, leading to the formation of 
a ^^C pocket. This mixing is needed in order to reproduce the 
enhancements observed in AGB stars of heavy elements pro- 
duced by slow neutron captures (the ^ process), with neutr ons 
released by the '^C(a,n)'*0 reaction dGallino et al.lll998l) . In 
the '^C pocket, the Mg isotopic composition can be altered 
by neutron captures. However, for IM-AGB models, we pre- 
viously demonstrated jKarakas et al.[ l2006bl) that the presence 
of a '^C pocket does not change the Mg isotopic ratios at the 
stellar surface because a reasonable value for the mass of the 



pocket is too small, ^ 10 Mq, to produce any effect on the 
overall Mg nucleosynthesis. 

Other important uncertainties in the stellar modeling re- 
gard how convectio n and c onvective borders are treated. 
IVentura & D'Antonal ( l2005bl) use a "Full Spectrum of 
Turbulence" (FST) prescription for convective regions, while 
we use Mixing Length Theory (MLT) with a-l.l. The FST 
approach leads to a smaller efficiency of the TDU and higher 
HBB temperatures. We cannot yet make a direct comparison 
with these models as detailed nucleosynthesis is only available 
at the rnqment for stars of Z = 0.001 with FST. As discussed by 
lHerwidll2004 gravity waves or convective extra mixing at the 
base of the envelope and at the boundaries of the thermal pulse 
can affect the amount of TDU m ass and also the temperature in 
the intershell. iHerwi ^ ('2004a'bl() employs a special scheme to 
allow for some time-dependent hydrodynamical overshoot at 
convective boundaries. He typically finds very deep TDU and 
also that the convective thermal pulse penetrates inside the C-O 
degenerate core. Models are available for Z = 0.0001, well be- 
low the metallicities we are considering here. Also in this case 
it remains to be analysed in detail if these models could provide 
a match to the composition of grain OC2. 

Finally, we conclude our discussion of model uncertain- 
ties by noting that rotation and magnetic fields are not in- 
cluded in our models. These typically affect the occurrence 
of extra-mixing phenome na and have been studied for stars 
of masses lower (see e.^. | penissenkov & VandenBergl 12003 : 
iTalon & Charbonnell 1200,^)7 or higher ("see e.g. iHeger et all 
200(1 iMaeder & Mevnell l2005t lYoon & LangeR l2005l) flian 
the IM star we are considering here. To our knowledge, how- 
ever, there are no published models of massive AGB stars with 
rotation. Rotation and magnetic fields are not ID phenomena 
and simplifications and parameterizations are required to input 
these physics into ID stellar structure codes. These simplifi- 
cations will result in extra parameters and uncertainties above 
those already included in the simulations. For these reasons 
it is also valuable to compute models without rotation, even 
though the effects of these phenomena on the evolution of mas- 
sive AGB stars should be studied in the future. 

3.3. The case for a low-mass low-metallicity AGB star 

A LM-AGB stai- of mass ^ 2.5 to 3.5 Mq and low metallic- 
ity, Z ^ 0.004 to ^ 0.008 can also produce the ^^Mg excess 
shown by OC2 because the He intershell temperature and the 
TDU efficiency are high enough in these stars. As an exam- 
ple, in Fig. jS] we present different models for a 2.5 Mq star 
of Z=0.004. As in the other figures, the predictions presented 
in this plot are calculated by including the abundance of ^^Al 
multiplied by the chosen factor of 25, however in these mod- 
els the contribution at mass 26 comes mostly from ^''Mg. In 
two runs of the post-processing of this model we artificially in- 
cluded the mixing of protons into the He intershell at the end 
of each TDU episode and we tested two different values for the 
extent in mass of the PMZ: 0.001 and 0.002 Mq. More ^^Mg 
is produced when the pocket is included because of neutron 
captures. Since this star is just not hot enough for HBB to oc- 



8 



M. Lugaro et al.: Peculiar spinel grain 0C2 



O 
oi O 
-> O 



CD 

■a o 



"T 



T" 



M=2.5, Z=0.004, no PMZ 

M=2.5, Z=0.004, Mp„2 = 0.001 M^^ 

M=2.5, Z=0.004, Mp„2 = 0.002 M,^ 




1000 



delta (^''Mg+25x^''AI)/^^Mg 

Fig. 8. The Mg predictions are shown for a 2.5 Mq Z = 0.004 
model with different choices of the extension in mass of the 
partial mixing zone (PMZ). 



cur we have to assume that some extra-mixing process at the 
base of the convective en velope, such as the c ool bottom pro- 
cessing (CBP) studied bv lNoUett et all ( l2003l) . is at work dur- 
ing the AGB phase, so that enough ^^Al is produced and the 
observed 5(^^Mg/^'^Mg) is matched. In order to reproduce the 
5(26Mg/24]y[g-) Q(-2 a 26A1/27A1 ratio in the range 0.06 to 
0.036 is required, wh ich can be produce d with CBP tempera- 
tures ^ 54 million K dNollett et al.ll2003h The low i**0/"'0 ra- 
tio in OC2 is consistent with efficient CBP and the '^O/'^O ra- 
tio would be 0.001 15 at the CBP temperature derived from the 
26^^27 ratio. Detailed calculations for CBP at low metallic- 
ities are required to test our present analysis before any strong 
conclusion can be drawn. We plan to include this process in our 
nucleosynthesis code in the future. 

3.4. Jhe Fe and Cr isotopic composition 

In AGB stars the Fe isotopic composition is modified be- 
cause the heavy Fe isotopes are produced by neutron cap- 
tures occurring during the convective thermal pulses when the 
^^Ne(Q', «)^^Mg neutron source is activated. This production is 
rather i nsensitive to the neutr on captures occurring in the '^C 
pocket jLugaro et all l2004ah . The (5("Fe/^<'Fe) value, at the 
time when 5(25Mg/24Mg)^433, is ^ 80 in our IM-AGB models 
and ^ 370 in the 2.5 Mq Z - 0.004 model. There are a cou- 
ple of reasons for this difference: first, the integrated neutron 
flux increases with decreasing the metallicity as less material 
means that more free neutrons are available. Second, the dilu- 
tion factor of the intershell material carried to the envelope by 
the TDU is about a factor of two higher in IM- than in LM- 
AGB stars. Unfortunately, the large uncertainty of the mea- 
sured d(^^Fe/^^Fe) prevents us from determining which models 
represent the best match. 

We cannot directly compare the Cr isotopic composition as 
we do not have the Cr isotopes in the network used to com- 
pute the present models. We can draw some ge neral features 
on Cr is otopes in AGB stars considering Fig. 2 of lLugaro et alJ 
ll2004ah . where results are presented for LM-AGB models of 



solar metallicity. The ^(^^Cr/^-Cr) and 5("Cr/"Cr) are barely 
modified by neutron fluxes in AGB stars while the 5(^"*Cr/^^Cr) 
values are slightly higher than the 5(^^Fe/^^Fe) values for a 
given stellar mod el. Since we are in the process of extend- 
ing our network jKarakas et all l2006al) . using preliminary 
models we have checked that these results also apply to our 
IM-AGB models. In fact, for the 5 Mq Z = 0.02 model 
we obtain no changes in (5(^°Cr/^^Cr) and (5(^^Cr/^^Cr), while 
6C'^CTr-CT)^nO, compared to 5("Fe/5^Fe)=71. For the 6.5 
Mo Z = 0.012 model we obtain 6C'^Ct/^^Ct)^40, at pulse num- 
ber 25. Our IM-AGB models would match the ^(^'^Cr/^^Q-) 
value measured in OC2, while the 2.5 Mq Z - 0.004 model 
would produce 6(^'^Cr/^~Ci) values close to the 2cr upper limit 
ofOC2. 

Finally, we note that in low metallicity stars the initial 
6{'Ci/^'Ci) could be largely negative because of the effect of 
the Galact ic Chemical Evolut ion of the Cr isotopes (see dis- 
cussion in IZinner et all Eo05l) . This represents a hint against 
a LM-AGB star of low metallicity as the parent star for OC2, 
even though detailed models for the evolution of the Cr isotopes 
in the Galaxy should be performed to confirm this analysis. 

4. Conclusions 

We have shown that the peculiar isotopic composition of O, Mg 
and Al in presolar spinel grain OC2 could be the signature of 
an AGB star of intermediate mass and metallicity close to so- 
lar (roughly higher than 0.008) suffering TDU and HBB, or of 
an AGB star of low mass and low metallicity (roughly lower 
than 0.008) suff^ering TDU and very efficient CBP While HBB 
occurs in our IM-AGB models, we do not simulate any extra- 
mixing in the LM-AGB model. Thus, the LM-AGB origin for 
grain OC2 has still to be tested against detailed models includ- 
ing CBP in low metallicity AGB stars. The large uncertainty in 
the Fe isotopic composition of OC2 does not allow us to de- 
termine which model better represents the parent star of OC2, 
but the Cr isotopic composition favors an origin in an IM-AGB 
star of metallicity close to solar. In this case, the model con- 
ditions to reproduce the composition of OC2 are well defined: 
a TDU mass > 0.05 Mq and/or a maximum He-intershell tem- 
perature > 360 million degrees, and a temperature at the base 
of the convective envelope in the range ^ 80 to 85 million de- 
grees. These conditions are satisfied by our 5 Mq, Z - 0.008 
and 6.5 Mq, Z = 0.02 models. Within this solution, we pre- 
dict that the '*0(/?,y)'^F and the '^©(p, a')'"*N reaction rates 
should be close to their lower and upper limits, respectively. It 
remains to be seen if the proposed rate changes are consistent 
with other constraints, such as future nuclear expriments, new 
grains similar to OC2, or the production of '^O in other stellar 
types. 

We note that using Salpeter's initial mass function, stars 
with mass from 5 to 7 Mq represent ^ 5% of all stars. Thus, one 
might expect such a proportion of presolar grains from IM- rel- 
atively to LM-AGB stars. In principle, any grain with "*0/'^0 
< 10""* and -*'Mg excesses is a candidate for an IM-AGB origin. 
There are 69 oxide grains in our database found by measure- 
ments of both '^O/'^O and "^O/'^O (thus representing an unbi- 
ased sample) and showing ^^Mg excesses. Of these, three have 
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i8q^i6q jQyygj- (jj^jj 10"^. Thus, at most ^ 4% of presolar oxides 
appe ar to have even the possib ility to have formed in IM-AGB 
stars. iBernatowicz et al.l(l2005h discuss detailed calculations of 
graphite and TiC grain growth in carbon stars. The ir Fig. 8 in- 
dicates that larger grains form in lower mass stars. iNuth et alJ 
propose a mechanism for the growth of crystals in AGE 
stars and red giants and show that only in low-mass (< 3 Mq) 
stars it is possible to form the large crystals found in primitive 
meteorites. Most of the oxide grains with both O and Mg mea- 
surements are larger than 1 micron, thus it is quite possible that 
our dataset is biased towards grains from low mass stars. 

Further analysis of presolar spinel grains may identify ad- 
ditional grains with isotopic compositions similar to 0C2 and 
give more precise measurements of their Fe and Cr isotopic 
compositions. This will provide the opportunity to test the find- 
ings of the present work as well as the possibility to constrain 
massive AGB models using presolar grains. 
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